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ABSTRACT 



Aims. The aim of this paper is to study the basic equations of the chemical evolution of galaxies with gas flows. In particular, we 
focus on models in which the outflow is differential, namely in which the heavy elements (or some of the heavy elements) can leave 
the parent galaxy more easily than other chemical species such as H and He. 

Methods. We study the chemical evolution of galaxies in the framework of simple models, namely we make simplifying assumptions 
about the lifetimes of stars and the mixing of freshly produced metals. This allows us to solve analytically the equations for the 
evolution of gas masses and metallicities. In particular, we find new analytical solutions for various cases in which the effects of 
winds and infall are taken into account. 

Results. Differential galactic winds, namely winds carrying out preferentially metals, have the effect of reducing the global metallicity 
of a galaxy, with the amount of reduction increasing with the ejection efficiency of the metals. Abundance ratios are predicted to 
remain constant throughout the whole evolution of the galaxy, even in the presence of differential winds. One way to change them 
is by assuming differential winds with different ejection efficiencies for different elements. However, simple models apply only to 
elements produced on short timescales, namely all by Type II SNe, and therefore large differences in the ejection efficiencies of 
different metals are unlikely. 

Conclusions. Variations in abundance ratios such as [0/Fe] in galaxies, without including the Fe production by Type la supemovae, 
can in principle be obtained by assuming an unlikely different efficiency in the loss of O relative to Fe from Type II supernovae. 
Therefore, we conclude that it is not realistic to ignore Type la supernovae and that the delayed production of some chemical elements 
relative to others (time-delay model) remains the most plausible explanation for the evolution of a-elements relative to Fe. 

Key words. Galaxies; abundances - Galaxies: evolution - Galaxies: ISM 



1. Introduction 

When a stellar system is formed out of gas, the dying stars 
enrich the remaining gas so that new stars are metal-enriched 
compared to their ancestors. Simple models of chemical evo- 
lution have been used since decades to interpret the evolution 
of the chemical composition of stars and gas. Pioneering works 
have determined the evolution of the metallicity in the Galaxy 
(Schmidt 1963, Searle & Sargent |T972l Tinslev [19741 Pagel & 
Patchett 1975 ) by making very simplifying hypotheses about the 
initial mass function (IMF), the lifetime of stars and the mixing 
of chemical elements with the surrounding interstellar medium 
(ISM). 

These models remain useful guides for understanding the 
main chemical properties of galaxies, but it was soon realized 
that they cannot give a complete picture of galaxy evolution. In 
particular, simple closed box models cannot reproduce the distri- 
bution of metallicity of long-living stars in the galaxy (G-dwarf 
problem, Lvnden-Bell |1975[ Rana |199"T| and references therein). 
The most obvious solution to the G-dwarf problem is to relax the 
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hypothesis that the solar neighborhood evolved as a closed box 
and to allow for infall of gas (Tinsley 1980). Increasingly com- 
plex models of chemical evolution including gas flows (infall, 
outflow, radial flows), but always assuming the instantaneous re- 
cycling approximation (IRAfl have been proposed, covering a 
large variety of galactic morphologies (Larson |1976t Clayton 
|19f8l Edmunds |T990l Koppen [1994 ; Martinelli & Matteucci 
I2OOOI 1. These models consider the global metallicity Z and ne- 
glect metallicity-enriched outflows, so they cannot be used to 
predict variations of abundance ratios. Pagel & Tautvaisiene 
(1998 ) have calculated abundance ratios using analytical mod- 
els, but they had to assume that chemical elements such as iron 
are ejected with a fixed time delay A after the time of star forma- 
tion. This is a rough approximation since it is known that the pro- 
genitors of SNela span a very wide lifetime (see e.g. Matteucci 
& Recchi 2001). Only detailed numerical models relaxing IRA 
and including the chemical enrichment from Type la SNe (e.g. 
Matteucci & Greggio 1986) allow one to follow in detail the evo- 
lution of single elements. Matteucci & Greggio d 19861 1 showed 
in detail the effect of the time-delay model, already suggested by 
Tinsley (1980) and Greggio & Renzini ( 119831 ). in particular the 
effect of a delayed Fe production by Type la SNe on abundance 
ratios involving a-elements (O, Mg, Si). 



' IRA assumes that stars with masses larger than 1 M© die instanta- 
neously whereas those below 1 M© live forever (Matteucci l200U 
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Outflows and galactic winds are produced when the thermal 
energy of a galaxy (produced by energetic events such as SN 
explosions and stellar winds) is enough to unbind a large frac- 
tion of the ISM (Mathews & Baker |1971"] l. The existence of dif- 
ferential winds, namely galactic winds in which the metals (or 
some of them) are ejected out of the parent galaxy more efli- 
ciently than the other elements, has been first hypothesized and 
applied to the chemical evolution of galaxies by Pilyugin ( 1993) 
and Marconi et al. ( |1994| l. In particular, Marconi et al. (1994) 
suggested that the chemical properties of gas-rich dwarf irreg- 
ular galaxies, most remarkably the spread in H/He vs. O/H and 
N/O vs. O/H, can be more easily explained by assuming differ- 
ential winds, in which products of SNell (mostly a-elements and 
one third of iron) can escape the potential well of a galaxy more 
easily than other elements formed in low and intermediate mass 
stars. Hereafter we will call selective winds the diff'erential winds 
in which different metals have diff'erent ejection efficiencies. The 
existence of differential and selective winds has been theoreti- 
caUy confirmed by several detailed chemodynamical simulations 
(Mac Low & Ferrara [T999l D'E rcole & Brighenti 1999; Recchi 
et al. 120011 Fujita et al. I2004l l and it has been observationally 
confirmed by determining the metallicity of outflows in starburst 
galaxi es (M artin et al. 120021 Ott et al. 2005). Recently, Salvadori 
et al. (I2OO8I 1 have suggested differential winds coupled with re- 
laxation of IRA in massive stars to explain the steep decline of 
the [O/Fe] ratio observed in the stars of Sculptor, without invok- 
ing any Fe production from Type la SNe. 

In this paper we apply the hypothesis of differential and se- 
lective winds in the framework of simple models of chemical 
evolution with IRA. The chemical evolution of galaxies is a com- 
plex phenomenon that can be properly treated only by means of 
detailed numerical codes, possibly taking into account also the 
hydrodynamical evolution of the system, as we have explained 
above. However, in spite of the oversimplifications introduced 
by simple chemical evolution models, the results they produce 
are very linear and simple to interpret. They can therefore illus- 
trate in a direct way the effect of differential winds on the chem- 
ical evolution of galaxies, as well as the effect of relaxing IRA. 
Moreover, for clarity purposes, we show also the results of de- 
tailed numerical models (Lanfranchi & Matteucci |2003t |2004| l 
as a comparison. 

The plan of the paper is as follows: in Sect. 2 we summa- 
rize the assumptions of the simple models of chemical evolution 
and the results of models including gas flows (infall and out- 
flow). In Sect. 3 we introduce the assumptions necessary to solve 
the equations of chemical evolution of galaxies with differential 
winds and produce some illustrative result. In Sect. 4 we study 
the abundance ratios of some chemical elements predicted by 
models with differential and selective winds. Finally, we discuss 
some implications of the results and draw the main conclusions 
in Sect. 5. 



2. Summary of simple model containing gas flows 

2.1. Standard solutions 

As it is well known, the so-called Simple Model of chemical evo- 
lution is based on the following assumptions: 

1. the system is one-zone and closed, namely there are no in- 
flows or outflows. 

2. The initial gas is primordial (no metals). 

3. The IMF is constant in time. 

4. The gas is well mixed at any time (instantaneous mixing). 



5. Stars more massive than 1 Mq die instantaneously; stars 
smaller than 1 Mq live forever [instantaneous recycling ap- 
proximation or IRA). 

These simplifying assumptions allow us to calculate analytically 
the chemical evolution of the galaxies. Once we have defined the 
fundamental quantities: 

Xoo 
(m - MR)(f>(m)dm, (1) 

(where (p(m) is the IMF and Mr is the mass of the remnant) and 



yz 



1 r°° 

Y^T^ J mpz,m<l>im)dm. 



(2) 



(where pz.m is the fraction of newly produced and ejected metals 
by a star of mass m), the well known solution of the closed box 
model can be easily found: 



Z = 3;zln(//-'), 



(3) 



where ji is the gas fraction MglM,ot. This result is obtained by 
assuming that the galaxy initially contains only gas and has 
the remarkable property that it does not depend on the partic- 
ular star formation history the galaxy experiences. To obtain eq. 
Q one has to assume that yz does not depend on metallicity. 
Nucleosynthesis can indeed depend on metallicity, especially for 
secondary elements. However, for O and Fe the metallicity effect 
is not so important (Matteucci l2001l Garnett l2002l l. therefore we 
can neglect it for the purposes of our work. In the general case 
in which yz - f(Z) the evolution of the metalhcity is governed 
by the following equation: 



dw 

7m 



ln(A^-'). 



(4) 



Relaxing the first of the assumptions of the simple model, 
we get the models including gas flows, also known as leaky box 
models. Analytical solutions of simple models of chemical evo- 
lution including infall or outflow are known since at least 30 
years (Pa gel & Patchett [19751 Hartwick 1976; Twarog [19801 
Edmunds 1 19901) . Here we follow the approach and the terminol- 
ogy of Tinsley (11980!) and Matteucci ( 20011) . namely we assume 
for simplicity linear flows (we assume gas flows proportional to 
the star formation rate (SFR)). Therefore, the outflow rate W{t) 
is given by: 

Wit) ^ A(l - R)m, (5) 

where iff{t) is the SFR, and the infall rate A(f) is given by: 

A(t) ^ A(l - R)il,(t). (6) 

Here A and A are two proportionality constants. The first as- 
sumption is justified by the fact that, the larger the SFR is, the 
more intense are the energetic events associated with it (in partic- 
ular supernova explosions and stellar winds), therefore the larger 
is the chance of having a large-scale outflow (see e.g. Silk 2003). 
A proportionality between A(t) and t//{t) can be instead justified 
by the fact that the infall of gas provides a continuous reservoir 
for the star formation. A careful analysis of this assumption is 
given in Sect. 13.41 
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Fig. 1. Evolution of the metallicity as a function of fi for the 
closed box model (solid line) and models with gas flows: A = 2 
and A = (dotted line); A = Q and A = 2 (short-dashed line); 
A - 3 and A = 1 (long-dashed line); A - I and A = 3 (dot- 
dashed line). The metallicity of the infalling gas Za is assumed 
to be 0. 



The system of equations we need to solve is therefore: 



dt 

dt 
dM. 



^{K-A){\-R)ilj{t) 
= (K-A-l){l-Rm) 



(7) 



^ £^ = (1 _ R)^i,(t)[hZA +yz-{A+ 1)Z] 



where Mz is the mass of metals (Mz - Z ■ Mg) and Za is the 
metallicity of the infalling gas. For the moment, Za is assumed 
to be a constant, although it could change with time due to the 
pollution of the intracluster medium (ICM) by means of the met- 
als ejected through the galactic winds. This aspect will be con- 
sidered in detail in Sect. 3.2 and Sect. 3.3. With some algebra, 
this set of equations can be worked out, yielding: 



AZa + yz 
A 



l_^(A_^)_(A_^_l)^-ij 



-1 A-,i-] 



(8) 



To find this solution we have made use of the initial conditions 
Z(0) = 0, Mg(Q) - M,o,{Q) - Mgo- It is worth pointing out 
that the explicit dependence of Z on time is hidden in the time 
dependence of fi (see also Sect. 13.4b . 

We show in Fig.[T]the evolution of Z/yz versus jj in the case 
Za = for different values of A and A. It is worth reminding 
that yu = 1 at the beginning of the evolution of the galaxy and its 
value decreases with time. Therefore, the time axis is inverted 
compared to the yu-axis in this figure and in similar ones. We 
can see that all the models with gas flows produce metallici- 
ties smaller than the ones predicted by the closed box model 
(solid line). This is in agreement with the first and the third of 
Edmund's theorems (Edmunds 1990), since we have assumed 
that the infall gas is pristine (Za - 0). Fig.|2]shows the effective 
yields y^jf - Zj ln(yU"') (normalized to yz) of the same models 



Fig. 2. Normalized effective yield as a function of yU. Notation 
as in Fig.[T] 

presented in Fig[T] y^/y is often used in literature to interpret data 
and to analyze to what extent the chemical evolution of a galaxy 
deviates from the behavior predicted by the closed box model 
(Garnett 2002', Tremonti et al. I2004I I. We can notice from eq. 
([8]l that the solution diverges for models in which /I > A when 
ju approaches 0, therefore, we show Z/yz and ye fflyz versus yu 
only until yii ~ 0.1. Indeed, Prantzos & Aubert (1 19951 1 pointed 
out that IRA cannot be applied in situations in which low yu are 
obtained, because the amount of gas ejected by low mass stars 
at late times can substantially modify the final abundances, also 
for elements produced on short timescales by massive stars. We 
notice here also that, in models in which the infall rate is larger 
than the outflow rate, not all the values of yU are allowed. Solving 
the first two equations of the system O, we get: 



A-i- 1 
^ K-A 



1 M, 



^,0 



K-AM„ 



(9) 



In the models in which K> A + \, Mg is always increasing (eq. 
|7]l, therefore yU ranges between 1 and a minimum value: 



A-^ 



(10) 



For models in which A + I > K > A there is no yUmin but there 
is a upper limit reachable by the gas mass which is given by 
Mg^Um - MgfiliA + 1 - A). For this reason we plotted in Figs. [1] 
and |2] the model with A - Q and A = 2 (short-dashed line) and 
the model with A - \ and A = 3 (dot-dashed line) only for 

yU > yUmin = 0.5. 



3. Simple models with differential winds 

3.1. The basic solution 

We present in this section the chemical evolution of galaxies in 
which a differential wind is assumed, namely in which the met- 
als are more easily channelled out of the parent galaxy than the 
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pristine gas. The easiest way to consider a differential wind in the 
framework of simple models of chemical evolution is to assume 
that the metallicity of the gas carried out in the galactic wind is 
proportional to the metallicity of the ISM with a proportional- 
ity constant larger than one. If we define Z" as the metallicity of 
the outflowing gas, this condition implies that Z" = aZ with the 
ejection efficiency a > 1. In the metallicity budget (third equa- 
tion in (|7]i) we have to assume that the negative term due to the 
galactic wind is given by 'W(t)Z° - aA(\ - R)\lj{t). A similar ap- 
proach has been used by Salvadori et al. (|2008) but in the frame- 
work of a semi-analytical model for the evolution of galaxies. 
With our simple approach we are able to determine analytical 
expressions for the evolution of Z, which allow us to understand 
more clearly the effect of galactic winds on the chemical evolu- 
tion of galaxies. 

The set of the equations we have to solve in this case is very 
similar to (|7]i, with the only difference given by the metallicity 
budget equation, which we modify as follows: 



dMz 
dt 



= (1 -R)m{^A+yz - {Aa + 1)Z]. 



The solution of this new set of equations is given by: 



Z = 



A^A +yz 

K + (a-\)A. 



r 

1 _[(A_^)_(A_^_ 1)^-1]- 



(11) 



(12) 



It is straightforward to see that we can obtain eq. dHJ in the case 
a = 1 (i.e. in the case in which the galactic wind is not dif- 
ferential). We remind here that this solution is valid only for pri- 
mary elements. The evolution of secondary elements (namely el- 
ements synthesized from the metals present in the stars at birth) 
in the framework of simple models with gas flows has been ex- 
tensively studied by Koppen & Edmunds (1 19991 1 and we refer 
the reader to this paper to find out more about it. 

Assuming again for simplicity an infall of pristine gas (i.e. 
Za - 0) we can see how the chemical evolution of the galaxy 
varies with the ejection efficiency a. In Figs. [3] and |4] we show 
Z/yz vs. yu for different values of a and gas flow parameters {A, 
A) taken as (3, 1) and (1, 3), respectively. The first couple of val- 
ues can be attributed to objects in which galactic winds have a 
prominent role. This is usually the case in dwarf galaxies, since 
the potential well of these objects is very shallow and large-scale 



outflows must be common (Dekel & Silk 1 19861 but see also 
Skillman 1 19971 1. The second set of parameters {A - \, K - 3) 
can more properly refer to spiral galaxies, although large out- 
flow episodes in massive spirals are not very likely. In general, 
the effect of differential winds can be more clearly distinguished 
in models in which the outflow is prominent compared to the in- 
fall. For this reason, we will pay from now on more attention on 
the models for which A > K, taking /I = 3, A = 1 as reference 
values. Our aim in this paper is not to describe a prototypical 
dwarf galaxy, therefore the combination of parameters (A, A) - 
(3, 1 ) is purely indicative and serves as a benchmark of our mod- 
els with differential winds. 

As expected, the differential wind determines a significant 
reduction of the metallicity of the galaxy. It can be shown ana- 
lytically that, for each value of a > 1, eq. ( fT2l l produces metal- 
licities systematically smaller that the ones obtained from eq. 
([8]). Moreover, it is evident from these plots the asymptotic trend 
of Zjyz as // approaches or ju„,„,. From eq. ( fT2b it is easy to 
see that, in both cases, Zjyz asymptotically tends to the value 
[(a - \)A + A]"' provided that Za - 0, therefore Zjyz decreases 
almost linearly with a when ji is small and o- » 1 . 
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Fig. 3. Evolution of the metallicity as a function of ji for simple 
models with differential winds, /i = 3, A = 1 and Za - 0. Plotted 
are models with the differential wind parameter a = 2 (dotted 
line), 5 (short-dashed line), 10 (long-dashed line) and 50 (dot- 
dashed line). For reference, also the model with a - 1 (normal 
wind) is plotted (solid line). 
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Fig. 4. As in Fig.[3]but for /i = 1 and A = 3 



3.2. Instantaneous ICM mixing models 

In our calculations we have so far assumed that the metallicity 
of the infalling gas Za is constant in time. This assumption is not 
necessarily true, in particular when differential winds are taken 
into account. In fact, the galactic winds can pollute the surround- 
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ing medium and the infalling gas can have a metallicity increas- 
ing with time. The increase of the iron abundance between red- 
shift z ~ 1 and z ~ in the ICM has been recently suggested 
(Balestra et al. I2007I I and can be accounted for in numerical 
models in which the chemical evolution of the ICM is affected by 
galactic winds and stripping (Calura et al. I2007I I. A simple hy- 
pothesis we can make about the variation of the ICM metallicity 
with time is that the metals, expelled out of the galaxy through 
galactic winds, instantaneously and uniformly mix with the sur- 
rounding ICM. This approach has been used by detailed chem- 
ical evolution models of the ICM (e.g. Matteucci & Vettolani 
11988 ; Moretti et al. I2003I I. Of course the instantaneous mixing 
of the outflowing gas with the whole ICM is rather unlikely. 
Nevertheless, this assumption is useful for illustrative purposes. 
This hypothesis leads to a new set of equations in which, be- 
sides Za, we have to consider also the variation, although small, 
of the ICM mass Mjcm due to galactic winds and infall. These 4 
equations are: 



^ = (1 - Rm)[AZA +yz- [Aa + 1)Z] 

dM; 



dt 

M,. 
dt 



(13) 



^^^^ = (1 -R)i//it)(A-A), 



where Mz^ is the mass of metals in the ICM (Mz^ - Za ■ Mjcm)- 
In order to solve this set of equations we have to define the ra- 
tio between gas mass and ICM mass - Mg/MicM and a new 
parameter k - Mg^/MicMiO), where M/cm(0) is the initial ICM 
mass. Once we have defined the values Xg and Xh such that: 



X-' = 



and 



A-^- 1 



(;t H- 1)(^ - A) H- 1 (k+\)iA-A) + l' 



XI 



k{A - A) 



k(A-A- 1) 

-1 . 



{k+l)iA-A) + l (yt+ 1)(^- A) + 1 
the solution of the system of equations (fTSl l is given by: 



yz 



AaXt' -(A-AXa- 1)L 



\-x. 



(14) 



(15) 



(16) 



It is straightforward to see that eqs (fT6] l and (fT2l) produce the 
same results if the infall is negligible (i.e. if A ~ 0). Considering 
again a model with A = 1 and /I = 3 we can show the trend of 
Zjyz in the case a - 5 for different values of k. We notice first 
that, unlike /i, rj does not monotonically decrease, but its varia- 
tion depends on the values of A and A. In the present example, 
since A > A, Mg decreases with time and Mjcm increases with 
time, therefore rj ranges between k at the beginning and during 
the late stages of the evolution of the galaxy. We show in Fig.|5] 
the evolution of Z/yz as a function of ri/k for k -0.1 (solid line) 
and k - 0.5 (dotted line). The model with k = 0.5 predicts larger 
abundances because the metals carried out by the galactic wind 
are diluted with less ICM, therefore Z^ is larger. These results 
cannot be compared to the results of eq. ( fT2] i in a straightforward 
manner since they are expressed in different units. However, we 
show in this plot the asymptotic trend of the simple differential 
wind model (short-dashed line in Fig. |3]l for jj approaching 
(dashed line). This comparison allows us to notice that, as ex- 
pected, eq. JTSl i predicts larger metallicities compared to models 
with Za = 0. However, as long as A; < 1, the diff'erences are small 
(of the order of ~ 0. 1 - 0.2 dex). 



0.08 



0.06 - 



N0^04 
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0,2 
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0.6 

rj/k 



0.8 



Fig. 5. Chemical evolution of models in which the metallicity 
of the ICM varies with time (instantaneous ICM mixing models) 
and with initial mass ratio k (see text) equal to 0. 1 (solid line) 
and 0.5 (dotted line). Also shown in this plot is the asymptotic 
trend of the simple differential wind model for // approaching 
(dashed line). For all models, A - 3, A - 1 and c = 5. 



We can also calculate the ICM metallicity, which turns out 
to be: 



(17) 



In the case we are considering more closely, A = 3 and A = 1 , 
therefore the ICM mass increases with time and Xg ranges be- 
tween 1 and (l+k). The larger value Z^/Z can attain is therefore: 



Za 
Z 



1 -(l +/t)"- 
which is smaller than 1 provided that: 



a-1 



(18) 



(19) 



Once the ratio Za/Z between ICM and ISM metallicity is obser- 
vationally known, eq. (fTsT l can be used to constrain a. 

3.3. Galactic fountain models 

A special case of variable infall metallicity is represented by the 
situation in which Za - Z" = aZ, namely the metallicity of the 
infalling gas is set to be always equal to the one of the galac- 
tic wind. This condition implies therefore that the very same gas 
that has been driven out of the galaxy by energetic events can 
subsequently rain back to the galaxy, due to the pull of its grav- 
itational potential. This kind of duty cycle is well known in as- 
trophysics and it has been named galactic fountain (Shapiro & 
Field 1 19761 Bregman |1980l ). In order to solve the chemical evo- 
lution in this special case, we have necessarily to assume that 
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0.1 - 



0.05 




Fig. 6. Comparison of the chemical evolution of a galactic foun- 
tain model (solid line) with a model with infall of pristine gas 
(Zyi = 0, dotted Une). For both models, /I = 3, A = 1 and a - 5. 

A > A since the reservoir for the infall gas is given by the gas 
expelled out of the galaxy through galactic winds. The solution 
of the chemical evolution in the presence of galactic fountains is 
given by: 



yz 



(A - A)ia - 1) 



1_[(A_^)_(A_^_ 1)^-1] 



U-A)(a-l) 
A-A- 1 



■(20) 



A comparison of this solution with the solution of eq. (fTZb 
with infall of pristine gas (Za = 0) is shown in Fig. |6] For 
both models plotted in this figure we have assumed A = 3, 
A = 1 and a = 5. As expected, the galactic fountain model 
predicts larger vales of Z. These values exceed also the metallic- 
ities attained by the models with instantaneous mixing of the 
ICM (Sect. 3.3) shown in Fig. |5] This result is expected be- 
cause the galactic fountain model maximizes the enrichment in 
metallicity of the infalling gas, whereas in the models shown 
in Sect. 3.3 the metals escaping from the galaxy are diluted in 
a large amount of ICM. From eqs. ( fTSl l and ( |20] l we can also 
notice that the metallicity ratio tends asymptotically to the value 
[(a-l)A+A]/[{A-A)(a-l)] which is always larger than one (for 
the special choice of parameters it is equal to 13/8) but it tends 
to 1 for Of » 1 . In our example, only an infall of pre-enriched 
gas, whose metallicity is larger than 5/8 yz can therefore pro- 
duce a metallicity in the galaxy larger than the one attained by 
the galactic fountain model. From eq. (l20l i we can also see that, 
in the special case in which A = A, Z/yz tends to the solution of 
the closed box model (eq. |3ll. This is due to the fact that, in the 
framework of the simple models of chemical evolution, having 
outflow and infall with the same rate and metallicity or not hav- 
ing gas flows at all is formally the same. Indeed, the fountains 
take a finite and non-negligible time to orbit around and fall back 
to the galaxy (Spitoni et al. [20081 ). This implies a delay in the 
mixing of metals in the ISM, which conflicts with the fourth as- 
sumption of the simple models of chemical evolution (Sect. 2), 
therefore only detailed numerical models can ascertain the effect 



of this delay on the chemical evolution of galaxies (Spitoni et al. 
2008, in preparation). 

3.4. Limits of fhe approximations of iinear fiows 

In the previous sections we have assumed that both outflow and 
infall rates are proportional to the SFR and we have seen that in 
this way we obtain solutions in which ifr(t) cancels out. In this 
section we analyze how reliable are these assumptions and what 
kind of results we obtain if we assume general infall and outflow 
laws (see also Edmunds 1 19901 Koppen & Edmunds 1999). First 
of all we notice that a correlation between infall rate and SFR 
arises naturally because the larger is the amount of infalling gas, 
the higher is the reservoir of gas inside the galaxy available to 
form stars. An example is given by detailed models of the chem- 
ical evolution of the Milky Way (Chiappini et al. 119971 Cescutti 
et al. 120071 Colavitti et al. 2008 ) in which two main phases of 
gas infall turn out to produce two distinct episodes of star forma- 
tion (although the similar behavior between the two rates does 
not imply that their ratio is constant). We will show in this sec- 
tion that, under reasonable assumptions, a constant ratio between 
SFR and infall rate can be attained during the late evolution of 
the considered galaxies. On the other hand, a proportionality be- 
tween the outflow rate and the SFR is quite realistic and it has 
been shown both observationally (e.g. Heckman 120021) and the- 
oreticafly (e.g. Silk 2003 ). 

Therefore, let us assume for simplicity that eq. Q still holds 
and that the infalling gas is pristine (e.g. = 0). The extension 
of our calculations to the case in which W{t) is a generic function 
is straightforward. If we assume a generic infall law A(t), the 
system of equations we need to solve is: 



^=A(r). 



: A(t) ■ 



dt 
dMf. 

IF 

dM; 



-A(l-Rm) 

(A+m-R)m 



(21) 



['I^^(l-R)^(t)[yz-(aA+l)Zit)l 



subject to the usual initial conditions Z(0) = 0, Mg(0) - 
MioiiO) - Mgfi. In this case, however, it is not possible to solve 
this system of equations canceling out i//(t). Therefore, at vari- 
ance with what we have done so far, we must assume some 
dependence of i//(t) on Mg (Schmidt law). We assume a linear 
Schmidt law (e.g. i/'(f) = SMg(t)), because it is the only formu- 
lation for which the results can be expressed analytically. For 
Schmidt laws with different exponents detailed numerical mod- 
els are required. The equation for the evolution of the gas mass 
can be simply solved, yielding: 



Mg(t) = e 



-(A+l)(l-R}Sl 



A{u)e 



(i+l)(l-R)S« 



du + M„o 



(22) 



Substituting this expression into the third equation of the system 
(ISTT i also the evolution of the metallicity as a function of time 
can be calculated. The result is the following: 



Z(f) = (1 -R)Syz[lit)] ' I(u)du, 



where 



/(f) = e 



rl 4,„iJ/i+l)(l-fi)5H 



(23) 



(24) 



We consider now as an example an exponential infall law, 
study what kind of metallicity evolution it produces and compare 
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the results with those presented in Sect. |3.U bearing in mind that, 
in the present case, we have to define the value of some constants 
(Mg o, S as well as the constants defining the infall rate), whereas 
in the previous case it was enough to define the proportionality 
constant A. We have: 



Ait) = Ke-'"' = Ke 



(25) 



where f, is a characteristic infall timescale and t = f/f,. We can 
calculate K assuming that the initial gas mass o is propor- 
tional to the total amount of infalling gas through a proportion- 
ality constant y (namely we assume that Mg o - 7 Jo Mt)dt). In 
this way, we obtain K - Eq. ( l22l i can be therefore solved: 



MJt) 



+ Kie 



-(A+l){l-R)StiT 



(26) 



where Ki is a constant (see the appendix for its exact value). 

We can now evaluate the ratio A(f) — ,, "1^,, , in order to test 
if and under which conditions we can still assume a proportion- 
ality between the SFR and the infall rate A(t). In Fig.Qwe show 
the resulting A(f) for diff'erent values of y (0.2, 1 and 5) and f, (1, 
1.5 and 5 Gyr). In order to draw these curves we have assumed 
/I = 3 (in agreement with what assumed in Sect. 13. Il l and 5 = 1 
Gyr"'. This is an upper limit for the star formation efliciency 
in dwarf galaxies (see. e.g. Lanfranchi & Matteucci I2003t . R 
is assumed to be 0.26, in agreements with the nucleosynthetic 
calculations of Woosley & Weaver (1 19951 1 assuming a Salpeter 
IMF. A small value of y implies that the galaxy initially does 
not have a large amount of gas and the star formation is mostly 
triggered by the infalling gas, therefore the infall rate initially 
prevails over the SFR. A large f, implies a large infall rate at 
late times, therefore A(f) tends asymptotically to a large value. 
In general the behavior of A(f) depends on a complex interplay 
between y, f,, S and A (see also Koppen & Edmunds 1999) but 
we can see that under reasonable assumptions the models show a 
first phase (lasting a few Gyr) in which A(f) varies with time, but 
then it stabilizes to a constant value. As a general result, only if 
we assume infall timescales much shorter than the star formation 
timescale (i.e. the inverse of S) A(t) shows significant variations 
over the whole evolution of the galaxy, although these variations 
are milder if a large value of A is assumed. For instance, if we 
take 5 = 0.1 Gyr"', a constant A{t) is attained by assuming ei- 
ther large values of A or large values of f, or values of y such 
that Ki ~ 0. If we assume a quadratic Schmidt law, the infall 
timescales required to obtain a constant A(t) tend to be larger 
than in the case of a linear Schmidt law. We can conclude that a 
linear flow can be a reasonable approximation of the late evolu- 
tion of a galaxy provided that the infall timescale is of the same 
order of magnitude of the star formation timescale. 

Now we can answer the following question: does an early 
phase of variable A(f) significantly affect the evolution of the 
metallicity? To do that we have to solve eq. ( l23T l. The resulting 
evolution of Z/yz is given by: 



Z(t) _K2+ K3e^" - [Ki + Ki)e-^^-^ 



(see the Appendix for the expressions of the involved constants). 

We show in Fig. [8] the evolution of Z(t)lyz as a function of 
time for different values of y. In these models we assume f, = 1 .5 
Gyr because this value leads to a A close to 1 (see Fig. |7]l. To 
compare these models with the results shown in Sect. 13.11 we 




5 10 

t (Gyr) 

Fig. 7. Ratio between the infall rate and the SFR as a function 
of time for models with an exponential infall and different f, and 
y (see text). Heavy lines: f, = 5 Gyr; normal lines: f, = 1.5 Gyr; 
light lines: f, = 1 Gyr Solid lines: y = 0.2; dotted lines: y = 1; 
dashed lines: y = 5. 

combine the first two equations of the system (|7|l to get yu(f) in a 
model with linear flows and a linear Schmidt law. By means of 
eq. ( fT2l i we then obtain the evolution of Z as a function of time, 
which turns out to be: 



Z(f) 1 - [(a - ^ + l) - (a - ^)e(^+i-A)(i-«)5'j A 

yz 



A-K»-ru 



A + {a-\)A 



(28) 



We plotted also this curve in Fig.[8]for A = 1 (long-dashed line). 
The difference with the models with exponential infall laws is 
very small and mostly due to the fact that the asymptotic value of 
A(f) for the exponential infall models is not exactly 1 . Therefore 
we can conclude that the results we have shown in Sect. 13. li do 
not change substantially if we assume an exponential infall in- 
stead of an infall rate proportional to the SFR. This demonstrates 
also that the major source of error in the solutions of the simple 
models is the IRA assumption rather than the assumption of lin- 
ear flows. 

4. The abundance ratios 

From eq. ( fT2l i it is straightforward to determine the abundance 
ratio between two chemical elements / and j assuming that the 
yields do not depend on metallicity (see also Sect.|2]l, namely: 



(27) - 



Z; Ji + AZa 



AZ, 



(29) 



This abundance ratio is therefore determined uniquely by con- 
stants which depend on the IMF and on the properties of the 
infalling gas. In the models with variable infall metallicities, 
Zi/Zj is simply given by the yields ratio yt/yj (see eqs. [T6land 
|20] i. With our simplifying assumptions, the abundance ratio must 
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Fig. 8. Evolution of the metallicity as a function of time for 
models with an exponential infall and different y. Solid line: y = 
0.2; dotted line: y = 1; short-dashed line: y = 5. For all the 
models it is assumed f,- = 1.5 Gyr. Also plotted (long-dashed 
line) the evolution of a model in which a linear infall (with A - 
1) is assumed. For all models a - 5. 



remain constant throughout the whole evolution of the galaxy. 
Also the abundance ratios of the ICM (see eq. [TtI i must re- 
main constant with (and therefore with time). The constancy of 
abundance ratios in the framework of simple models of chemi- 
cal evolution was also pointed out by Matteucci & Chiappini 
(I2005I I. in which the authors showed that variations of the abun- 
dance ratios can be accounted for only by relaxing the IRA and 
by properly taking into account stellar lifetimes. Here we outline 
that this result holds also if we consider differential winds. 

By neglecting the contribution of SNela, it is expected that 
the [O/Fe] ratio does not change with time (and therefore with 
[Fe/H]) since, in this case, both elements are produced by SNell, 
on similar timescales. Indeed, detailed models of the evolution 
of the Galaxy (Chiappini et al. '1997) show a mild decrease of 
the predicted (and observed) [O/Fe] as a function of [Fe/H] even 
without considering SNela, due to the variation of the yields with 
stellar mass, but this can account for only ~ 0.2 dex of the [O/Fe] 
variation. A much more significant reduction of the [O/Fe] can 
be attained by taking into consideration SNela, which provide 
most of the Fe of the ISM and explode with a delay compared to 
the SNell. This model (also called time-delay model) remains the 
most plausible explanation of the variation of [O/Fe] vs. [Fe/H] 
in various galaxies (Matteucci & Greggio II986I Lanfranchi & 
Matteucci 120031 1. In particular, Lanfranchi & Matteucci ( 120031 1 
interpreted the steep decline of the [a/Fe] ratios at low [Fe/H] (at 
[Fe/H] — 1.6 dex), observed in dwarf spheroidals of the Local 
Group (e.g. Shetrone et al. 1200 11 Venn et al. |2004| l as due to the 
slow SFR acting in these systems, coupled with very efficient 
galactic winds. In this situation, in fact, the Fe production from 
Type la SNe occurs when the [Fe/H] is still low. At this point, the 
galactic winds extract both a-elements and Fe from the galaxy 
but while O is less and less produced because of the decrease 



in the star formation due to the wind, Fe continues to be pro- 
duced by Type la SNe. Unfortunately, SNela cannot be treated 
by means of simple model (SNela arise from CO white dwarfs 
in binary stars, with non-negligible lifetimes, therefore the IRA 
cannot be applied). Because of that, it is not possible to solve an- 
alytically chemical evolution equations which include SNela and 
this is the reason why most people make use of numerical sim- 
ulations. However, IRA provides a good approximation of the 
evolution of a elements, promptly released by SNell. Recently, 
Savadori et al. (2008) proposed an alternative scenario in which 
the drop in the [O/Fe] is due to strong differential winds, com- 
bined with the relaxation of the IRA, without invoking SNela. 
In their model, the star formation is dominated by an early, very 
short (~ 100 Myr) burst. However, a non-negligible fraction 
of SNela seems to appear also within 0. 1 Gyr in galaxies (see 
Mannucci et al. 2005 ; 2006), therefore neglecting SNela in their 
model could be incorrect. This might be the reason why these 
authors do not reproduce accurately the high metallicity tail of 
the metallicity distribution function of Sculptor (namely the stars 
with [Fe/H] 1; see their figs. 4 and 6). On the other hand, ne- 
glecting Type la SNe appears consistent with the stellar lifetimes 
they have assumed. These lifetimes predict 150 Myr for a 8 M© 
star, instead of 30^0 Myr as commonly adopted. These long 
lifetimes, together with the short duration of the main episode 
of SF and the relaxation of IRA help in producing the decrease 
of [O/Fe] since in this case the evolution of abundance ratios re- 
sembles that of an instantaneous burst in which the [O/Fe] at a 
given time t follows the ratio of O and Fe yields of the star dying 
at that time. 

One way to attain a variable [O/Fe] ratio in the framework 
of simple models of chemical evolution is by assuming selective 
winds, namely that the two elements have different ejection effi- 
ciencies. It has been suggested since a long time that a-elements 
can be channelled out of a galaxy more easily than other ele- 
ments (Marconi et al. 1994; Martin et al. 2002) but Recchi et al. 
(.2001) showed that, during the very first burst of star formation 
experienced by a dwarf galaxy, iron can be expelled more easily 
than oxygen. The hypothesis of different ejection efficiencies a 
for O and Fe is therefore a reasonable assumption that can be 
tested by our models. 

We assume two different values for these ejection efficiencies 
ag and Ofe and, by means of eq. ([T2l i. we calculate how [O/Fe] 
varies with [Fe/H]. The results would not change significantly 
if we assume eqs. ( fT6l ) or eq. ( |20l ) instead of eq. ( fT2b . In Fig. |9] 
we show for instance the comparisons of models for which (ag, 
Ufe) - (5, 2) (solid line), (2, 5) (dotted line) and (2, 2) (dashed 
line). We show the results for our reference set of parameters 
(A - 3 and A = 1) in the left panel. To study the dependence of 
these curves on the inflow and outflow parameters, we also plot 
the results for A = 1 and A = 3 (right panel). In order to calculate 
the yields, we have assumed a Salpeter IMF between 0.1 and 40 
Mq and the nucleosynthetic prescriptions of Woosley & Weaver 
([1995) for solar metallicity (case B). For the solar abundances 
we adopt the Anders & Grevesse 1 1989 J values. We use this set 
of solar abundances in order to be consistent with the choice of 
Woosley & Weaver's ( 119951 ) yields. 

As expected, as ffie galaxy ages (and, therefore, [Fe/H] in- 
creases), [O/Fe] bends down if the oxygen ejection efficiency is 
larger (solid line), whereas it starts increasing if af^ > ag (dot- 
ted line). The model with A = 1 and A = 3 shows the same 
behavior but with reduced variations of the abundance ratios. 
This is due to the fact that the difference between the metallic- 
ities attained for a - 2 and a = 5 in the case (A, A) = (3, 1) 
is larger than in the case (A, A) = (1, 3) (see Figs. |3] and |4]i. Of 
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Fig. 9. [O/Fe] vs. [Fe/H] for models characterized by different 
ejection efficiencies for O and Fe ap and Ofe . (ao, ape) are (5, 
2) (solid line), (2, 5) (dotted line) and (2, 2) (dashed line). {A, A) 
are (3, 1) (left panel) and (1,3) (right panel). 



course, the model in which the ejection efficiencies are the same 
(dashed line) does not show any variation of [O/Fe] with [Fe/H] 
(see eq.|29]l. The [O/Fe] attained by this model (~ 0.2) is lower 
than the [O/Fe] of halo stars and of the most metal-poor stars of 
Local Group dwarf galaxies, whose chemical enrichment should 
have been dominated by SNell ejecta. Indeed, it is known since 
quite some time that Woosley & Weaver's yields overpredict iron 
(Timmes et al. |1995t Chiappini et al. 1997). As we have said, 
the model with A = 3 and A = 1 can be representative of dwarf 
spheroidal galaxies. Indeed, in this model the metallicity remains 
significantly lower than the solar value and, if O is ejected more 
easily than Fe, a knee at [Fe/H] — 1.5 appears, in agreement 
with observations (see e.g. Venn et al. |2004| l. We point out once 
more that we are neglecting SNela therefore, in the framework 
of simple models, O and Fe are both produced by massive stars, 
on similar timescales. As a consequence, they should share a 
similar fate and have comparable ejection efficiencies. Recchi et 
al. (12004 ). by means of detailed chemodynamical simulations, 
found similar values of ap and a/r^ when complex star forma- 
tion histories are assumed. Therefore, simple models predict that 
variations in the [O/Fe] vs. [Fe/H] plot cannot be attributed to 
differential winds and SNela must necessarily be considered in 
order to properly interpret this diagram. 

Would this conclusion still hold if we relax IRA and take 
properly into account stellar lifetimes? In order to test that, we 
have performed a series of numerical simulations of the chemi- 
cal evolution of a dwarf spheroidal galaxy whose structural pa- 
rameters resemble those of Local Group such as Carina, Draco, 
Ursa Minor, Sculptor. In particular we have run a simulation 
very similar to the "standard" model of Lanfranchi & Matteucci 
( |2003i l. The standard model is adjusted to reproduce the average 
[or/Fe] and neutron capture abundance ratios of six local dwarf 
spheroidals. It is assumed that the galaxy forms through the fast 
and continuous infall of pristine gas until a mass of ~ 10'^ is 
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Fig. 10. [O/Fe] vs. [Fe/H] for numerical models of a stan- 
dard dwarf spheroidal galaxy. Solid line: wind model like in 
Lanfranchi & Matteucci (2007); dot-short-dashed line: the same 
model without SNela; long-dashed line: the same model with- 
out SNela and without galactic winds; short-dashed line: model 
without SNela in which the ejection efficiency of ff-elements is 
enhanced. Data are from Bonifacio et al. (|2000i 2004); Shetrone 
et al. (|200T[ [2 003]l; V enn et al. |20 04l Sa dakane et al. (2004); 
Fulbright et al. (120041 1; Geisler et al. (I2005I I. For clarity purposes, 
the dot-short-dashed line has been shifted of 0.0 1 dex since oth- 
erwise it would have been overlapped with the long-dashed line. 



reached. The star formation consists of one long episode (8 Gyr) 
with low efficiency (v = 0.05 Gyr"') and is affected by very in- 
tense galactic winds (with a rate 2 times higher than the SFR). 
As the galactic wind develops, it removes a large gas fraction 
thus decreasing the SFR. This effect, coupled with the low star 
formation rate and the injection of Fe into the ISM by SNe la, 
gives rise to the knee observed in [(r/Fe] and [r-process/Fe]. The 
models we show now contain non-selective winds. We plotted 
the [O/Fe] vs. [Fe/H] resulting from the standard model in Fig. 
[TO](solid line). We have then run the same model without SNela 
(dot-short-dashed line) and without SNela and without galactic 
winds (long-dashed line). We have also considered a model with- 
out SNela in which the ejection efficiency of a-elements is en- 
hanced (selective winds) (short-dashed line). Similarly to what 
we have seen in Fig.|9] if we neglect SNela, there is a negligible 
variation of the [O/Fe] ratio as a function of [Fe/H] (dot-short- 
dashed and long-dashed lines) and a knee can be attained only 
in the presence of a-enhanced winds, namely winds in which O 
is ejected out of the galaxy more efficiently than Fe (ap = 2a f^; 
short-dashed line). The turnoff point of this curve corresponds to 
the moment in which the galaxy experiences the galactic wind. 
This plot shows that only a model which takes properly into ac- 
count SNe la and galactic winds is able to produce the knee in 
[ff/Fe] observed in several dwarf spheroidal galaxies. 
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5. Discussion and conclusions 

In this paper we have studied general properties of simple chem- 
ical evolution models of galaxies with infall and differential 
winds, namely with a selective loss of metals from the parent 
galaxy. We have seen that, assuming gas flows proportional to 
the SFR, the equations of the evolution of gas mass, total mass 
and metallicity Z of gas can be solved analytically, leading to a 
simple expression of the variation of Z as a function of the gas 
mass fraction ji and of the ejection efficiency of metals a (which 
is defined as the ratio between the metallicity of the galactic 
wind and the metallicity of the ISM). As expected, the result- 
ing metallicity of the galaxy turns out to be systematically lower 
than the metallicity attained by models with non-differential 
winds. In particular, the metallicity of models with differential 
winds tends asymptotically (for low values of fi, namely dur- 
ing the late evolution of the galaxy) to a constant value which 
decreases almost linearly with a. We have checked that, for ex- 
ponential infall laws and star formation rates proportional to the 
gas mass (linear Schmidt law) a constancy between the star for- 
mation rate and the infall rate arises quite naturally during the 
late evolution of the galaxies, provided that the star formation 
timescale does not differ significantly from the infall timescale. 
In the first few Gyr the ratio between infall rate and star for- 
mation rate can vary quite significantly but, in spite of that, the 
resulting evolution of the metallicity as a function of time does 
not change appreciably compared to the case with linear flows. 

We have also analyzed models in which the metallicity of the 
infalling gas varies with time. A realistic case of infalling gas 
with increasing metallicity can be envisaged in the interaction 
galactic winds/ICM. Indeed, this interaction must certainly lead 
to some pollution of the gas external to the galaxy, therefore to 
a constant increase of the ICM metallicity, in particular when 
differential winds are considered. In order to take it into account, 
we have made two extreme assumptions: 

- the metals carried out of the galaxy instantaneously and uni- 
formly mix with the surrounding ICM {instantaneous ICM 
mixing). 

- The metallicity of the infalling gas is set to be equal to the 
metallicity of the outflow, representing a physical situation 
in which the very same gas carried out of the galaxy by out- 
flows (or a fraction of it) falls back to the galaxy due to the 
pull of its gravitational potential {galactic fountain). 

Real physical situations should lie between these two extreme 
assumptions since the distribution of metals in the ICM is ex- 
pected to be patchy (Domainko et al. 2006; Sanders & Fabian 
12006 ; Tornatore et al. 2007). The metallicities predicted by in- 
stantaneous ICM mixing models are of course larger than mod- 
els in which the infalling gas has a pristine chemical composi- 
tion. The magnitude of this difference depends on the relative 
amount of ISM and ICM masses. Under normal circumstances, 
in which the ICM mass is larger than the ISM one, the differ- 
ence in the final metallicity is only ~ 0. 1 - 0.2 dex. A larger ISM 
metallicity is obtained by galactic fountains models. This result 
due to the fact that galactic fountains maximize the chemical en- 
richment of the infalling gas. This metallicity tends however to 
the one attained by the model with Za = if very large metal 
ejection efficiencies (a » 1) are assumed. 

We have then studied the metallicity ratios produced by dif- 
ferential wind models. We have seen that the models with vari- 
able infall metallicities (both instantaneous ICM mixing and 
galactic fountain) simply predict that abundance ratios between 



chemical elements are given by ffie ratio of their respective stel- 
lar yields (if these do not change with metallicity). Assuming 
a fixed infall metallicity Za Q this ratio also depends on the 
value of Za and on the infall parameter A, but nevertheless it 
does not change during the evolution of the galaxy. However, 
we stress here that simple models can be properly applied only 
to SNell products, whose release is very prompt. Iron is for in- 
stance mostly produced by SNela, whose explosion timescales 
range from a few tens of Myr to several Gyr, depending on the 
star formation history of the considered galaxy (Matteucci & 
Recchi I2001I I. therefore only the fraction of iron produced by 
SNell (~ 1/3 of the total) can be treated by the simple models. A 
constant O/Fe ratio, as predicted by simple models is therefore 
expected because, in this case, both chemical elements are pro- 
duced by massive stars, on similar timescales and they should 
therefore have a similar evolution. However, O/Fe ratios are not 
constant in galaxies and larger ratios are observed in older, more 
metal-poor stars. The most natural explanation of this trend is 
that the delayed production of iron by SNela reduces the initially 
high O/Fe ratio of the oldest stars {time-delay model; Matteucci 
& Greggio [T986] l. 

Simple models can produce variable abundance ratios either 
if we assume different ejection efficiencies for different chem- 
ical species (selective winds) or an IMF variable in time. Both 
assumptions should be quite "ad hoc". We have tested the case of 
selective winds for the O/Fe ratio showing that, if the O ejection 
efficiency is significantly larger than the Fe one, O/Fe decreases 
with time (and with Fe/H). We have confirmed this result also by 
means of detailed numerical simulations in which IRA is relaxed 
showing that, neglecting SNela, the only way to get a knee in 
the [O/Fe] vs. [Fe/H] relation is by assuming a-enhanced galac- 
tic winds. In a model suitable for dwarf spheroidal galaxies, in 
which the outflow rate exceeds the infall rate, the O/Fe is almost 
constant up to [Fe/H] — 1.5, then it bends down. This behav- 
ior is consistent with the observation of many dwarf spheroidal 
galaxies of the Local Group (Venn et al. 2004 and references 
therein). If instead iron is ejected more efficiently than oxygen, 
the opposite happens, with O/Fe ratios increasing with time. On 
the other hand, an iron ejection efficiency larger than the oxy- 
gen seems plausible, as it has been shown by chemodynamical 
simulations of dwarf galaxies experiencing violent bursts of star 
formation (Recchi et al. I2001tl2002l l. However, also these mod- 
els predict O/Fe ratios decreasing with time (and with [Fe/H]) 
because SNela were not neglected and their contribution is more 
important than selective winds in determining the evolution of 
the O/Fe ratio. Moreover, Recchi et al.'s model found larger iron 
ejection efficiencies precisely because of the delayed explosion 
of SNela. These SNe explode in a medium already heated and 
diluted by the previous activity of SNell, therefore SNela prod- 
ucts can be very easily channeled along the galactic funnel. If 
we neglect SNela, O and Fe are produced by the same stars and 
they should therefore share the same fate. 

Our main conclusions can be summarized as follows: 

- The differential winds have the effect of reducing the global 
metallicity of a galaxy. The magnitude of this reduction in- 
creases with the ejection efficiency of metals. 

- This remains true also if we allow a variable metallicity of 
the infalling gas. Larger galactic metallicities are obtained if 
we assume that the infall has the metallicity of the outflowing 
gas (galactic fountain model). 

- Simple models with constant (i.e. non metallicity- 
dependent) yields predict constant abundance ratios, 
even in the presence of differential winds or variable 
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infall metallicities. The same occurs for detailed numerical 
models. 

- One way to get variable abundance ratios in simple models 
is by assuming selective winds, namely differential winds 
with different ejection efficiencies for each chemical ele- 
ment. However, simple models apply only to chemical ele- 
ments promptly released on short timescales, therefore large 
element-to-element variations of the ejection efficiencies are 
not very likely. Abundance ratios predicted by simple mod- 
els of chemical evolution are then likely to remain constant 
during the evolution of the galaxy. 
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6. Appendix 

Here we give the explicit value of some constants we have used 
in Sect. 13^ 



El = iA+l)(l -R)Sti-l. 
E2 = A{a-l){l-R)Sti. 
Ki ^yEi-l. 



K3 



A(a- I)' 

(l-R)Sti 
(l-R)(l+aA)Sti-l' 



(30) 
(31) 
(32) 

(33) 
(34) 



